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CRITICAL COMPRESSIVE STRESS FOR OUTSTANDING FLANGES 

By Eugene B. LtjSdqtjist and E&bbidge Z. Stowell 



SUMMARY 

A chart is presented for the values of the coefficient in 
file formula for the critical compressive stress at which 
buckling may be expected to occur in outstanding flanges. 
These flanges are flat rectangular plates supported along 
the loaded edges, supported and elasticaUy restrained along 
one unloaded edge, and free along the other unloaded edge. 

The mathematical derivations of the formulas required 
for the construction of the chart are given. 

INTRODUCTION 

In the design of s tressed-skin structures for aircraft 
as well as in the design of compression members, it is 
desirable to know the compressive stress at which 
buckling occurs. In practice the structure is usually 
so imperfect or so eccentrically loaded that lateral 
deflection starts with the beginning "of loading. When 
lateral deflection starts with the beginning of loading, 
however, there is usually a very, pronounced increase in 
deflection at the critical compressive stress for which 
buckling would have occurred had the structure been 
perfectly straight and centrally loaded. The evalua- 
tion of this critical compressive stress for a flat plate, 
with certain conditions of edge support, is discussed 
in this paper. 

When a flat plate is loaded in compression, the two 
loaded edges are either simply supported or restrained in 
some manner. If the two unloaded edges 'are not sup- 
ported, the plate is considered to be a column. When 
one, or both, unloaded edges are also supported or 
restrained in some manner, the critical compressive 
stress is greatly increased over that for the plate as a 
column. That the compressive stress is increased when 
one, or both, edges are supported or restrained in some 
manner has been recognized for years. Because of the 
importance of the edge conditions, formulas based on 
the assumption that each edge of the plate is free, 
simply supported, or fixed have been employed in 
design. (See the s umma ry of these formulas given in 
reference 1.) . 

A study of the theory and the more reliable test 
data on the buckling of plate elements in s tressed-skin 
structures and compression members revealed the neces- 



sity for a more careful consideration of the edge condi- 
tions of plates than has been previously attempted. 
Accordingly studies were made of the critical compres- 
sive stress for I-, Z-, channel, and rectangular-tube 
sections in which proper consideration was given to the 
interaction between the individual parts of the cross 
section. (See references 2, 3, and 4.) In order to 
make the results of the work more generally applicable, 
studies were also made of the basic plate elements that 
comprise these sections. All the design charts resulting 
from this investigation were made available in 1938. 
The combination of the present paper with references 
2, 3, 4, and 5 is a more complete presentation of all this 
material. 

The basic element treated in this paper is a plate 
simply supported along the loaded edges, supported and 
elastically restrained against rotation alongone unloaded 
edge and free along the remaining unloaded edge. This 
basic element is representative of the outstanding flange 
on the I-, Z-, and channel-section columns. In reference 
5 is treated the basic element representative of the webs 
of these sections with elastic restraint along both 
unloaded edges. 

The mathematical derivations required for the investi- 
gation of the present paper are given in appendixes A 
and B. The results of practical use are given in the 
body of the paper. 

EVALUATION OF CRITICAL STRESS 

Within the elastic range. — Within the elastic range 
in which the effective modulus of elasticity is Young's 
modulus, the critical compressive stress fa- for a thin 
flat rectangular plate is expressed as (reference 6, p. 331, 
equation (214)) 

where 

k nondimensional coefficient that depends upon condi- 
tions of edge restraint and shape of plate 
E Young's modulus 
t thickness of plate 
p. Poisson's ratio 
b width of plate 

III 
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Beyond the elastic range. — When the pJata is 
stressed in compression beyond the elastic range, the 
effective modulus of elasticity for the plate is less than 
Young's modulus. If a single, over-all effective plate 
modulus riE is substituted for Young's modulus E, the 
critical stress, when the material of the plate is loaded 
beyond the elastic range, can be obtained from equa- 
tion (1). The nondimensional coefficient rj has a value 
that lies between zero and unity and is determined by 
the stress. For stresses within the elastic range, i?=l. 
For a more complete discussion and definition of tf, see 
reference 2. 

TltjEis substituted for Sin equation (1), the resulting 
equation cannot be directly solved for j tt . If the 
equation is divided by i\, however, f„/ri is given directly 
by the geometrical dimensions of the plate, Young's 
modulus E, and Poisson's ratio p. Thus 

>„_ k*>EP 
v 12(1 -flb* 



(2) 



For a given material, the relationship between/ er and 
fer/ri tends to be fixed by the compressive stress-strain 
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Fiouhi 1.— Variation of/«- with /«/<> for 24S-T aluminum alloy of minimum required 
properties. (When kipa/sq In., it- 1 and /«,-/,/,,) 



curve. This relationship is discussed in reference 2, 
where it is shown how probable relationships between 
j„ and/„/ij are obtained from the column, curve of the 
material because column curves are more readily avail- 
able than compressive stress-strain curves. The ques- 
tion is, therefore, what column formula should be used? 
Equations (8) and (9) of reference 7 define column curves 
that apply when the material just satisfies the mini- 
mum requirements of Navy Department Specification 
46A9a for 24S-T aluminum alloy. The relationships 
between / cr and/„/i; for this case are given in references 
2, 3, and 4 and in figure 1 of this paper. 

The 24S-T material delivered under specification 
46A9a almost always has properties that are better 
than the minimum required properties. The relation- 
ships between /„ and /„/ij for the average 24S-T 
material delivered are given in figure 2. This figure 
has been prepared in the manner described in reference 
2, the column curves for average 24S-T material as 
given in reference 8 being used. 



Figures similar to 1 and 2 of this paper may bo pre- 
pared for any material. The engineer using this paper 
must therefore decide whether the computation should 
be based on minimum required material properties or 
average material properties. 

Regardless of whether figure 1 or 2 is used, if the 
restraint against the rotation of the flange at its base is 
near zero and \fb is greater than approximately 2.5, 

2 



it is recommended that the curve ij 
For all other values of the restraint, 



- be used, 
the curve 



v= T ' should be satisfactory. In figures 1 and 2 the 

different equations mvolving r merely identify different 
curves that result from the relationships indicated. 
The value of r is E/E, the ratio of the effective column 
modulus for bending failure at the stress f er to Young's 
modulus. 

When the restraint against the rotation approaches 
zero, the i ? = T "^V^ curve is recommended in recognition 
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FiQUH2 S.— Variation of /„ with {J* for 248-T aluminum alloy of average propertl**. 
(When/^,<lfl.7klpa/sqln.,,-lanu/„-/«/i.) 

of the fact that the resistance of the plate elements to 
buckling arises largely from their torsional rigidity. 
The two curves recommended to show the relationship 
between j„ and/ C r/i? should be used until future experi- 
mental data indicate that different curves should bo 
used, 

EVALUATION OF * 

, The value of j er fa at which buckling occurs is given 
by equation (2), in which all of the quantities are known 
except the value of the coefficient k. The values of 
k can be obtained from figure 3; figure 3(b) is a portion 
of figure 3(a) plotted to a larger scale. In this chart, 
k is plotted against the ratio of the half wave length to 
the width X/6 for different values of a parameter e, 
termed the "restraint" coefficient. (In reference 9 
Trayer and March refer to « as the "fixity" coefficient. 
In this paper e is called the restraint coefficient to avoid 
confusion with the fixity coefficient c for columns.) 
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Fiau«*3.— General design chartering values of k. — — j^fizpjji 



The restraint coefficient c depends upon the relative 
stiffness of the plate and the restraining element along 
the side edge of the plate. The simplest conception 
of e is obtained when the restraining element, or stiff- 
ness, is assumed to be replaced by an elastic medium 
in which rotation at one point does not influence rota- 
tion at another point. For this type of restraining 
medium along the edge of the plate, 



l M & T 

within the elastic range, e=— ^- 

beyond the elastic range, e ^^=T 
where 



(3) 
(4) 



So stiffness per unit length of elastic restraining medium 
or moment required to rotate a unit length of 
elastic medium through one-fourth radian 

D flexural rigidity of plate, per unit length J_ i2(i— M «)J 

^ coefficient to allow for a decrease in D due to the 
■application of stresses beyond the elastic range 
Inasmuch as ij is a function of stress, its value for 
24S-T material can be obtained from figure 4 or 5, 
depending upon whether minimum required properties 
or average properties are being used. The values of 
t u T 2) and Tj also given in figures 4 and 5 occur in 
appendix A. 
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If S 0 is zero, e is also zero and the condition of simple 
support, or zero restraint, is obtained. If So is infinite, 
e is also infinite and the condition of a fixed edge or of 
infinite restraint is obtained. Therefore a variation of 
e from zero to infinity will cover all possible conditions 
of restraint at" the side edge of the plate. 

Figure 3 shows that for each value of e there is a 
value of X/& for which it is a minimum. Strictly, a 
whole number m of half wave lengths X must exist in the 
length of the platen. Hence 



Thus, to read a value of k from figure 3, it is necessary 
to substitute m=l, 2, 3, etc. in equation (5) until a 
value for \/b is obtained that gives the smallest value 
of it in figure 3. This smallest value of k is the one 
to be used in equations (1) or (2). This general pro- 
cedure will always give the correct value of k for use in 
equations (1) or (2) regardless of whether or not S 0 , 
and hence «, is a function of the half wave length X. 

For the special case in which S a , and hence t, is inde- 
pendent of the half wave length X, the general procedure 
described for obtaining a value for k can be used to 
construct a new chart, with the abscissa X/6 replaced 
by-alb. This new chart is given in figure 8. 

When S 0 , and hence e, varies with X or X/6, figure 6 
should not be used, but the general procedure as 
applied to figure 3 should be used to obtain the correct 
value of k for equations (1) and (2). 

EVALUATION OF S„ AND « 

Before it is possible to determine k from figure 
3 or 6, it is necessary first to evaluate the restraint 



coefficient e. The value of S 0 to be substituted in 
equation (3) or (4) will depend upon the characteristics 
of the structural member or members that provide tho 
restraint. In this paper it is assumed that the restraint 
is provided by a specially defined elastic restraining 
medium. As a result of this assumption, it has been 
possible to derive the general chart of figure 3, which 
is independent of the structure that provides tho 
restraint. 

The basic property of the elastic-restraining medium 
is that rotation at one point of the medium docs not 
affect rotation at another point of the medium. In 
many practical problems the elastic restraint is pro- 
vided by a stiffener, a plate, or some other structure for 
which rotation at one point affects rotation at another 
point. Consequently, the evaluation of S 0 in any 
given problem must take into account the effect of this 
interaction within the elastic restraining structure. 

The formula for So to be used in any given problem 
will depend upon the type of structural member that 
provides the restraint. Because this entire subject of 
restraint supplied to the side edge of a plate has been 
rather superficially treated in the literature, it is being 
made the subject of a series of papers by tho NACA, 
the first of which is reference JO. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
I Langley Field, Va., March 14, 1941. 



APPENDIX A 



SOLUTION BY DIFFERENTIAL EQUATION 



The procedure for obtaining the critical stress of a 
plate uniformly compressed along two opposite, simply 
supported edges is given in reference 6 (p. 337). In 
this method, which was also used by Dunn in reference 
11, the critical stress is found by solving the differential 
equation expressing the equilibrium of the buckled 
plate. The same method is applied in this paper to the 
case in which an elastic restraint against rotation is 
present along one unloaded edge of the plate while the 
other unloaded edge remains free to deflect and to 
rotate. For generality, the elastic restraint is assumed 
to arise from an elastic medium distributed along the 
unloaded edge; this medium has the basic property that 
rotation at one point within it does not influence the 
rotation at any other point. 

Figure 7 shows the coordinate system and the plate 
dimensions. The differential ^equation for the equilib- 
rium of a plate element is 



CA-1) 



where 

/ uniformly distributed compressive stress 

t thickness of plate 

w deflection normal to plate 

x longitudinal coordinate in direction of applied stress 

D flexural rigidity of plate, per unit length 

y transverse coordinate across width of plate 

ti, 1*, and t 3 coefficients equal to or less than unity 

In equation (A-l) the term ^(d'w/dx 2 ) is concerned 
with the external forces on the plate that cause buck- 
ling; whereas the term 



is concerned with the internal resistance of the plate to 
buckling. The terms involving r t and t, in equation 
(A-l) are concerned with the longitudinal and the 
transverse bending, respectively; whereas the term 
involving r t is concerned principally with the torsional 
stiffness. The coefficients t u r t , and r s allow for the 
change in the magnitude of the various terras as the 
plate is stressed beyond the elastic range. In the 
elastic range, n=T t =ri=l. 

The loaded edges are simply supported and are not 
displaced in the direction to. Of the several forms 




Figure 7. — Outstanding flange under edge compression. 

of the general solution of equation (A-l) the following 
form was selected as appropriate for this problem: 

w=((7 l cosh ^+ d sinh^-f- (k cos f+ C t sin fXos y 

(A-2) 

where 



and 



i2(i-jg)6y 



(A-5) 



Equation (A-2) satisfies the boundary conditions at 
the loaded edges and gives real values for both a and jS 
near the buckling stress /=/„. 

117 



118 



REPORT NO. 734 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



The values of the coefficients &, &, C t , and C 4 are 
to be found from the boundary conditions along the side 
edges of the plate. The value of X, the half wave length 
of the buckle pattern, is found from the condition that 
there must be an integral number of half wave lengths in 
the length a of the plate; thus 



m 



(A-6) 



where m=l, 2, 3, etc. 

In the elastic range, where T 1 =r 1 =r 3 =l, the values 
of a and 0 are 



(A-7) 



(A-8) 



The solution given by equation (A-2) was selected to 
satisfy the boundary conditions of no deflection and 
simple support (no moment) along the loaded edges. 
The boundary conditions along the unloaded side edges 
have also to be satisfied. The boundary conditions 
along the unloaded side edges are: 

(w),„ 0 =0 



(A-9) 
(A-10) 

(A-ll) 

(A-12) 



where So is the stiffness per unit length of the elastic 
restraining medium or the moment required to rotate a 
unit length of the medium through one-fourth radian. 



From equations (A-9) and (A-10) are obtained 



(A-13) 
(A-14) 



where 

€= 4g6 , (A _ lg) 

From equations (A-ll) and (A-12) are obtained 
Cjjjp sinh a + a t°^0i (P cosn «+2 cos 0) J 

-C^g sin p--J±_(p cosh a+g cos 0)J=O (A-16) 
cj^qa cosh a-f-^jqp0i(tfa sinh a—pfi sin 0)J 

— C^pP cos 0--j£l^( 2a sinh a-pp sin 0)J=O (A-17) 

p^-^J (A-18) 



where 



(A-19) 



The buckled form of equilibrium of the plate is 
obtained when the determinant formed by the coeffi- 
cients of Ct and C 4 in equations(A-16) and (A-17) equals 
zero. 

Thus, 

(a*+ 0*) (p*0 sinh a cos 0— g*a cosh a sin 0) 



+ «[(?>*+ 2*) «0 cosh a COS 0+2pga0 
+ (p'0*— gV) sinh a sin 0]=O 



(A-20) 



This equation establishes the critical compressive stress 
for an outstanding flange elastically restrained against 
rotation at one unloaded side edge. Thus equation 
(A-20) was used to establish the exact values of lc 
given in table I. 

The condition of simple support (no restraint) along 
the supported edge is described by f=0. For this 
special case, the problem is to find the smallest value of 
k?*0 that will satisfy equation (A-20) when e=0. A 
convenient method for determining this value of k is 
first to solve for e: 

_ (g*+0 l ) (p'0 sinh a cos 0-g'a cosh a sin 0) 
€ (p* +2 s )a^ cosh a cos 0 + 2p%aP+(p*p - g^a^sinh a sin0 

(A-21) 

(A-22) 
(A-23) 



When e=0, either 



or 



««+/3 2 =0 

p*P sinh a cos 0— <fa cosh a sin 0=0 



or 



(p*+q r )aP cosh a cos 0+2pga0 

+ (p*0*- g'a 2 ) sinh a sin 0= ~ (A-24) 

Equation (A-22) is true only if k=0, which can be 
true only if the compressive stress / is zero. Equation 
(A-24) applies only if k= <» , which can be true only if 
the compressive stress / is infinite. Consequently if a 
finite value of it: y^O forwhich e=0 exists, equation (A-23) 
must be satisfied. 

The special case of a fixed side edge (infinite restraint 
along the supported edge) is described by t= <» . 
Equation (A-21) shows that, if «= ~ , either 

a >+ i 8 , =o 0 (A-25) 

or 



p'0 sinh a cos 0— g*a cosh a sin f}=a> (A-26) 



or 



(P*+2 S )«0 cosh a cos 0+2j>ga0 

+ feV-aV) sinh « sin 0=0 (A-27) 

Equation (A-25) is true only for k= », which can be 
true only if the compressive stress / is infinite. Equa- 
tion (A— 26) cannot be true for a finite value of k. 
Hence if a finite value of k for which e= » exists, 
equation (A-27) must be satisfied. 



APPENDIX B 



SOLUTION BY ENERGY METHOD 



Because the exact solution of the differential equation 
given in appendix A does not lend itself to a direct 
calculation of k as in tlie case of the energy method of 
solution, an energy solution was made to aid in the 
construction of the chart of figure 3. The energy 
method gives approximate values for k, the accuracy of 
which depends upon how closely the assumed deflection 
surface describes the true deflection surface. 

The energy method as applied to the calculation of 
critical compressive stress is given in reference 6 (p. 327). 
The plate is stable when (\\+ Vi)~> Tand unstable when 
(T'i+FsXr, where J 1 is the work done by compressive 
forces on the plate, \\ is the strain energy in the plate, 
and V| is the strain energy in the elastic restraining 
medium along one side edge of the plate. The critical 
stress is obtained from the condition of neutral stability: 
T=\\+V, (B-l) 

If w is the deflection normal to the plate at any 
point x, y in the plane of the plate shown in figure 
7 and S 0 is the stiffness per unit length of the elastic 
restraining medium or moment required to rotate a unit 
length of elastic medium through one-fourth radian, 
then T, V u and \\ are given by the following equations 
(see reference 6, equations(199)and(201)and reference?), 
equation (73)): 



"di 1 by*. 



1 



dxdy (B-3) 



CM. 

In order to evaluate T, \\, and V s , it is necessary 
to assume a deflected surface w consistent with the 
boundary conditions. These boundary conditions at 



7 OS 893 0 -45-8 



the 1 side edges of the plate are, in the coordinate system 
of figure 7, 

(w) r . 0 =0 (B-5) 

When buckling occurs, a restraining moment will be 
applied to the plate along the edge y=0; the magnitude 
of the moment will depend upon the stiffness of the 
elastic restraining medium. If the elastic medium 
offers no restraint against rotation, this moment will be 
zero and the plate will swing about the edge y—0, as 
about a hinge. In this case the plate will remain essen- 
tially flat across its width. On the other hand, if the 
elastic medium offers infinite restraint against rotation, 
the plate will not rotate along the edge y=0 and the 
plate will deflect across its width into a shape similar 
to that for a cantilever beam. For any restraint of the 
elastic medium between zero and infinity the deflection 
curve across the width of the plate is taken as the sum 
of the straight line and the cantilever-deflection curve. 
In the direction of the length the usual sine curve 
indicated by the solution of the differential equation is 
used. Thus the deflection surface assumed for the 
plate is, in the coordinate system of figure 7, 

«=H+<(f)'+<f)+<l)Xf)1-? 

(B-9) 

where A and B are arbitrary deflection amplitudes and 
fli=— 4.963, o l =9 J 852, and ^=—9.778. These values 
of a u CE a , and a t were selected by taking the. proportion 
of two deflection curves that gave the lowest critical 
compressive stress for a fixed-edge flange for which 
/i=0.3. These two deflection curves were for a canti- 
lever beam with iateral uniform load and for a lateral 
load proportional to y. 
The condition 5=0 represents the case of a simply 
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supported or hinged edge at y=0. The case of A=0 
represents the condition of a clamped edge at y—0. 
The ratio ,A/B is therefore a measure of edge restraint 
and is related to the restraint coefficient e through the 
boundary condition given in equation (B-6). Sub- 
stitution of w as given by equation (B-9) in equation 
(B-6) gives 



where, by definition, 



4^6 



(B-10) 



(B-ll) 



Substitution of the value of B as given in expression 
(B-10) in the deflection equation (B-9) gives 

-^[(!)+<!)+<I)+<f)l C08 ? 

(B-12) 

Equation (B-12) shows how the shape of the deflection 
surface is affected by the restraint coefficient «. This 
equation is used in the evaluation of V u Vi, and T. 
Thus, 



\\=A- 



2b\ 



+3 



m 



(B-14) 



(B-15) 



where 



c 1 =-(l+a I +a,+o l ) ==0.79546 

c 3 =-(6+5a I +4a a +3a I ) =0.89395 

. ai J , 2a s . . 2gia 8 . a, 1 



c 4 = 



/ 



+^f^?+°f-')=0.04286 



144 



1 /400 . 



i l , +48a 1 + 120,'+ 16a I a 3 +4a,* 
+80a 1 +20a 1 +36aia 1 +12a,<z,^==0.56712 
1 /25 , 16 . . 30 . 9 . . 16 

0.17564 



+|ai s + 5a, + j a, + 4a 1 a J + 3^,)= ( 



CT=^i(5+9a 1 +8a 8 +7a t +4a l *+7a 1 a a +6a 1 a, 

+3a, , +5a 2 a,+2a 1 ')=0.19736 

c 8 =aaC,= — 2.3168 
Cj =a, 2 c 4 =4.0982 

It is permissible to substitute the values of T, T'i, 
and V t as given by equations (B-14) to (B-16) in 
equation (B-l) only when the applied stress / has its 
critical value/ er . After this substitution it is found that 

k**EP 



where 



1 , Cgt ,~c?~ 
3" r 2a," r 4a J t 



(B-16) 



(B-17) 



Equation (B-17) was used to calculate the values 
of k listed in the columns designated (a) of table I. 
With these values of £ as a guide, a number of correct 
values of k were obtained by satisfying equation (A-20) 
of appendix A. In this manner the errors in k as 
given by equation (B-17) were established at isolated 
points. From this knowledge of the errors, corrections 
were made to all the values of k given in columns (a) 
of table I. These corrected_values of .it, which are 
recommended, are listed in the columns designated (b) 
of table I. The recommended values of k were used 
in the construction of figures 3 and 6. 
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TABLE I. — VALUES OF k IN THE BUCKLING FORMULA FOR OUTSTANDING FLANGES ELASTICALLY 

RESTRAINED AT THE BASE 




•Values obtained from the energy method. 
• Recommended values. 

« Values obtained (com. the si act solution of the differential equation. 



